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Abstract. We give a precise mathematical formulation for the notions of a parametric 
patch and linear precision, and establish their elementary properties. We relate linear 
precision to the geometry of a particular linear projection, giving necessary (and quite 
restrictive) conditions for a patch to possess linear precision. A main focus is on linear 
precision for Krasauskas' toric patches, which we show is equivalent to a certain rational 
map on CP d being a birational isomorphism. Lastly, we establish the connection between 
linear precision for toric surface patches and maximum likelihood degree for discrete ex- 
ponential families in algebraic statistics, and show how iterative proportional fitting may 
be used to compute toric patches. 



Introduction 

Bezier curves and surfaces are the fundamental units in geometric modeling. There are 
two basic shapes for surfaces — triangular Bezier patches and rectangular tensor product 
patches. Multi-sided patches are needed for some applications, and there are several control 
point schemes for C°° multi-sided patches. These include the S'-patches of Loop and 
DeRose [17], Warren's hexagon [22], Karciauskas's M-patches [12], and the multi-sided 
toric Bezier patches of Krasauskas [14]. (Relationships between these and other patches 
are discussed in [13].) Parametric patches are general control point schemes for C°° patches 
whose shape is a polygon or polytope. They include the patch schemes just mentioned, as 
well as barycentric coordinates for polygons and polytopes [9j [211 123] . 

The success and widespread adoption of Bezier and tensor-product patches is due in 
part to their possessing many useful mathematical properties. Some, such as affine in- 
variance and the convex hull property, are built into their definitions and also hold for 
the more general parametric patches. Other properties, such as de Castlejau's algorithm 
for computing Bezier patches, come from the specific form of their Bernstein polynomial 
blending functions. Linear precision is the ability of a parametric patch to replicate linear 
functions. When the blending functions of a parametric patch correspond to the vertices of 
a polytope, these blending functions give barycentric coordinates precisely when the patch 
has linear precision. In this way, blending functions for a parametric patch having linear 
precision are barycentric coordinates for general control point schemes. 
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For us, a (parametric) patch is a collection of non-negative blending functions indexed 
by a finite set A of points in M. d , where the common domain of the blending functions is 
the convex hull A of A. A collection of control points in M indexed by A is used to define 
a map from A to M> e . The blending functions determine the internal structure of this map 
and the basic shape, A, of its image, while the control points determine how the image 
lies in M. e . Choosing the control points to be the points of A gives the tautological map, 
and the patch has linear precision when this tautological map is the identity on A. We 
show that every patch has a unique reparametrization (composing the blending functions 
with a homeomorphism of A) having linear precision (Theorem II. lip . This generalizes 
Theorem 8.5 of [IS], which was for toric patches. 

The blending functions of a patch may be arbitrary non-negative C°° functions. A patch 
is rational if it has a reparametrization having polynomial blending functions. For a rational 
patch, a map given by control points corresponds to a linear projection of a projective 
algebraic variety associated to its blending functions. Its unique reparametrization having 
linear precision has rational blending functions if and only if a certain canonical map 
defined on this variety is a birational isomorphism (Theorem 12 . 9[) . which implies that this 
variety has a maximally degenerate position with respect to a canonical linear subspace 
given by the set A (Theorem 12.101) . 

We apply this analysis to Krasauskas's toric patches, which are rational. The shape of 
a toric patch is a pair (^4, w) where A is a set of integer points in Z d and w is a collection 
of positive real numbers indexed by A. We show that iterative proportional fitting [5], 
a simple numerical algorithm from statistics, computes the blending functions that have 
linear precision. This algorithm was suggested to us by Bernd Sturmfels. It may form a 
basis for algorithms to manipulate these patches. 

Krasauskas [T5], Problem 3] asked whether any toric patches (besides the classical Bezier 
simploids [Hj) admit a rational reparametrization having linear precision. A toric patch of 
shape (A, w) corresponds to a homogeneous polynomial / = f^ w whose dehomogenization 
is the sum of monomials with exponents A and coefficients w. The toric patch admits a 
rational reparametrization having linear precision if and only if the toric differential 
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defines a birational isomorphism (Corollary 13. 13p . This analysis of linear precision is used 
in [11] to classify which toric surface patches can have linear precision. 

In Krasauskas's question about linear precision, he allowed the points A indexing the 
blending functions of toric patches to move within their (fixed) convex hull, keeping the 
same blending functions (See Example 13.161) . Our analysis of linear precision for general 
parametric patches will help to address that version of his question. 



In Section 1, we define parametric patches and show that every patch has a unique 
reparametrization that has linear precision. In Section 2 we show that if a rational patch 
has linear precision, then an algebraic variety we obtain from the blending functions has 
an exceptional position with respect to a certain linear space. We define toric patches 
in Section 3, and in Section 4 we explain how iterative proportional fitting computes the 
blending functions for toric patches that have linear precision. 
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1. Linear precision for parametric patches 

We interpret the standard definition of a mapping via control points and blending func- 
tions (see for example [131 §2]) in a general form convenient for our discussion. All functions 
here are smooth (C°°) where defined and real- valued unless otherwise stated. Let M> be 
the set of strictly positive real numbers and M> the set of non-negative real numbers. 

Let A be a finite set of points in M. d , which we shall use as geometrically meaningful 
indices. A control point scheme for parametric patches, or (parametric) patch for short, is 
a collection (3 = {/3 a | a £ A} of non- negative functions, called blending functions. The 
common domain of the blending functions is the convex hull A of A, which we call the 
domain polytope. We will always assume that A is full dimensional in that it has dimension 
d. We assume that the blending functions do not vanish simultaneously at any point of A. 
That is, the blending functions have no base points in their domain. 

A set {b a | a £ .4.} C M. e of control points gives a map F: A -> R e defined by 

(1.1) F(x) - E-eA&O*) 1 *- 



The denominator in (11.11) is positive on A and so the map F is well-defined. 

Remark 1.2. Positive weights, or shape parameters, {w a £ 1R> | a £ A} scaling the 
blending functions are often included in the definition of the map F. We have instead 
absorbed them into our notion of blending functions. See Remarks 11.71 and 13.31 @ 

The blending functions {/3 a | a £ .4.} are normalized if they form a partition of unity, 

= 1 - 

aeA 

As the denominator in (II. ip is (strictly) positive on A, we may divide each blending function 
by this denominator to obtain normalized blending functions. 

Redefine (3 a := ^ a . Then (3 a (x) > and Pa{%) = 1; 

for x £ A. For normalized blending functions, the formula (II. ip becomes 
(1.3) F(x) = ^/3 a (x)b a . 

ae-4 

We deduce two fundamental properties of parametric patches. 

• Convex hull property. For x £ A, F(x) is a convex combination of the control 
points. Thus F(A) lies in the convex hull of the control points. 

• Affine invariance. If A is an affine function on M £ , then A(F(x)) is the map with 
the same blending functions, but with control points {A(b a ) | a £ ^4}. 

A patch is non- degenerate if its blending functions are linearly independent. This implies 
in particular that F(A) cannot collapse into a point unless all the control points are equal. 
We make our key definition. 
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Definition 1.4. A patch {f3 a | a G A} has linear precision if for every affine function A(x) 
defined on M. d , 

A(x) = ^/5 a (x)A(a). 

That is, if the patch can replicate affine functions. This notion depends strongly on the 
positions of the points in A. ^ 

The tautological map r of a given patch is the map (11.31) when the control points are 
taken to be the corresponding points of A, 



(1.5) t(x) := £fl 



a (x) a. 



aeA 



By the convex hull property, t(A) C A. By affine invariance, the patch has linear precision 
if and only if r is the identity function, r(x) = x. We record this fact. 

Proposition 1.6. A patch has linear precision if and only if its tautological map is the 
identity map on the domain polytope A. 

Remark 1.7. This is a more restrictive notion of linear precision than is typically consid- 
ered. It is more common to define a patch to be a collection | i — 1, . . . , n} of blending 
functions, and this has linear precision if we have 



for some non-negative weights {wi \ i = 1, . . . ,n} and some points A = {a z | i — 1, . . . , n} 
whose convex hull is the domain polytope A of the blending functions. 

Our more restrictive definition of patches and linear precision, where we incorporate 
the weights and the points A into our definition of blending functions, allows us to derive 
precise criteria with which to study linear precision. Our intention is to employ these 
criteria to study patches in the generality in which other authors had worked. ^ 

Example 1.8. Fix an integer n > and let A :— {- | < i < n} be the set of n+1 
equally distributed points in the unit interval. The ith Bernstein polynomial is (3i := 
(^)af (1— x) n ~ l , which we associate to the point - £ A. These form a partition of unity, 

n n / \ 

i=0 i=0 ^ ' 

and are therefore normalized blending functions. Given control points bj G M 3 , for- 
mula (11.11) becomes 



i=0 

which is the classical formula for a Bezier curve of degree n in 
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This patch has linear precision. First, note that - ■ (™) = Then r(x) is 



i=0 



n 



i=l 



n-l 
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x^l-x)' 
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= x -J2 

3=0 

Example 1.9. Let A be the vertices of the hexagon, given below by its defining inequalities. 

l-y>0 



1 + x - y > 



1 + x > 




x > 



1 +|/ — x > 



1 + y > 

Let the blending function associated to a vertex be the product of the linear forms defining 
edges that do not contain that vertex. 

:= (l+x-y)(l-y)(l-x)(l+y-x) 

:= (l-y)(l-x)(l + y-x)(l+y) 

/3[0] := (l-x)(l + y-x)(l + y)(l + x) 

:= (l+y-x)(l + y)(l + x)(l + x-y) 

:= (l + y)(l + x)(l + x-y)(l-y) 

/3[_0] := (l+x)(l+x-?/)(l-y)(l-x) 

The normalized blending functions have denominator ip(x,y) := 2(3 + xy — x 2 —y 2 ), which 
is strictly positive on A. These normalized blending functions have linear precision. 



1.10) 



^/3 a (x,y)a = <p(x,y) 



In fact, these normalized blending functions are Wachspress's barycentric coordinates [21] 
for this hexagon. ^ 

When A is the set of vertices of A (so that A is in convex position), then normalized 
blending functions of a patch with linear precision are barycentric coordinates for A. By 
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this we mean in the sense of pi ETJ [23] : nonnegative functions indexed by the extreme 
points of A that have linear precision in the sense of Definition II .41 

A reparametrization of a patch {/3 a | a £ -4} by a homeomorphism if): A — > A is a 
new patch with blending functions (3 a o ip for a £ ^4. A patch is proper if its tautological 
map (11.51) is a homeomorphism. This condition is necessary for the patch to have linear 
precision. If a patch is proper, then reparametrizing it by the inverse r _1 of r gives a 
new patch whose tautological map is r o r" 1 , the identity function, and r _1 is the unique 
homeomorphism of A with this property. We record this straightforward, yet fundamental 
result about linear precision. 

Theorem 1.11. A proper patch has a unique reparametrization that has linear precision. 

Theorem II . 1 1 1 suggests that it will be fruitful to discuss patches up to reparametrization. 
In Section [21 we show that a natural geometric object X@ associated to a patch is invariant 
under reparametrizations and thus represents the patch up to reparametrization. 

A patch is rational if it has a reparametrization whose normalized blending functions 
are rational functions (quotients of polynomials). A patch has rational linear precision Hit 
is proper and its reparametrization having linear precision has blending functions that are 
rational functions. Such a patch is necessarily rational. We seek criteria that determine 
when a rational patch has rational linear precision. 

Theorem 1 1.111 and our preceding discussion concerns reparametrizations of a patch having 
linear precision. We may also alter a patch by multiplying its blending functions by positive 
real numbers, typically called weights, or by moving the points of A that are not vertices of 
the domain polytope A. This second change does not alter any map (11.11) given by control 
points, but it will change the tautological map r (ll.5p . Our analysis in the next section 
may be helpful in addressing whether it is possible to tune a given patch (using weights 
or moving the points in A) to obtain one that has rational linear precision. Example 13.161 
(due to Krasauskas and Karciauskas) shows how one toric patch may be tuned to achieve 
linear precision. 

2. The geometry of linear precision 

We introduce an algebraic-geometric formulation of patches to clarify the discussion in 
Section [1] and to provide tools with which to understand and apply Theorem 11.111 This 
leads us to discuss the relevance for geometric modeling of the subtle difference between 
rational varieties and unirational varieties. Lastly, we use this geometric formulation to 
give a geometric characterization of when a patch has rational linear precision. We recom- 
mend the text [3] for additional background on algebraic geometry. We first review linear 
projections, which are the geometric counterpart of control points. 

2.1. Linear projections. We consider M. e as a subset of the ^-dimensional real projective 
space MP* via the embedding 

R E 3 z i — ► [l,z] £ MP*. 
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A point [zo, Zi, ... , zp] G WP e lies in this copy of M. e if and only if zq ^ 0. In that case, the 
corresponding point of Mf is 

' Z\ z 2 Zi ' 



(2.1) 

w Zq Zq Zq 



Let A C R be a finite set of points and write MP for the real projective space whose 
coordinates are indexed by A. Let b = {b a G M £ | a G .4.} be a collection of control points, 
which we regard as points in MP £ so that b = {[1, b a ] G MP £ | a G A}. 

Given a point y = [y& | a G A] G MP" 4 , if the sum 



, a , G M m 



(2.2) ^y a -(l,b B 

is non-zero then it represents a point ir^y) in MF e . This map 7Tb is a linear projection 

(2.3) tt 6 : MP -4 ► R¥ e . 

We use a broken arrow > to indicate that ttj, is defined only on the complement of 

the set Ef, C MP" 4 where the sum (12.21) vanishes. This linear subspace is the center of 
projection. 

When the control points aflinely span WL e , the center of projection Ej, has codimension 
£ + 1 in MP" 4 . The inverse image 7r fe ~ 1 (x) of a point x G MP^ is the set H \ Eb, where 
H C MP^ is a linear subspace of codimension t containing Eb. If L is any ^-dimensional 
linear subspace of MP -4 that does not meet Eb, then 7^ maps L isomorphically to WP e . 
Identifying L with MP gives an explicit description of the map 7T b . If y G MP- 4 \ E b , then 



7r 6 (j/) = L n (j/, s 6 ) , 

the intersection with L of the linear span of y and the center Ej,- If y G L, then its inverse 
image under 7r fo is (y, E b ) \E b . 

For example, Figure [1] shows the effect of a linear projection 7r on a cubic space curve 
C . The center of projection is a line, E, which meets the curve in a point, B. 




Figure 1. A linear projection it with center E. 
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2.2. Geometric formulation of a patch. Consider ( 11. ip from the point of view of pro- 
jective geometry. Given a patch (3 = {f3 a (x) | a G .4} with A C M d , then 

(2.4) x i — ► \pjx) | a G .4] 

defines a map f3: A — ► MP- 4 whose ath coordinate is the blending function (3 a (x). 

The map /3 is unchanged if we multiply all blending functions by the same positive 
function <p(x) defined on A. Thus we obtain the same map if we use instead the normalized 
blending functions. The image Xp of A under (3 is a closed subset of MP -4 (but not Zariski 
closed!), as A is compact and the map (3 is continuous. It is non-degenerate (does not lie 
in a hyperplane) exactly when the patch is non-degenerate. 

Suppose now that we have control points b = {b a G Mr | a G .4} and consider the 
composition of the map f3 and the linear projection -k^ 

A X/3 C MP" 4 MP^. 

By our assumption on the positivity of the blending functions, formula (12. ip implies that 
the image lies in the standard copy of Mr, and there it is given by the formula (II. ip . Thus 
this image, F(A), is the image of Xp under the projection 7r& given by the control points. 

A reparametrization of the patch f3 by a homeomorphism ip : A — > A gives a different 
map with the same image in MP- 4 

?: A A x,. 

Thus X/3 is an invariant of the patch modulo reparametrization. 

If we take ^4 to be our set of control points where the element a G A is the control 
point associated to the blending function j3 a , then the resulting linear projection n A is the 

tautological projection, 

(2.5) ix A : MP" 4 ► MP d . 

Here, the coordinate points (e a | a G ^4) are indexed by elements of A and 7T4(e a ) = a G 
MP d . The tautological map r ( II. 5p is the composition 

A X p c MP^ ---^ MP d . 

The patch is proper if this map is a homeomorphism onto its image A C MP d , and the 
patch has linear precision if the composition is the identity map, by Proposition 11.61 

With these definitions, we have the following identification of the reparametrization 
having linear precision, which is a more precise version of Theorem 11.111 

Theorem 2.6. Suppose that Xp is a proper patch and ir^: Xp — > A is the tautological 
projection restricted to Xp. Then the blending functions for Xp that have linear precision 
are given by the coordinates of the inverse ofix^. 

Remark 2.7. This geometric perspective, where a patch is first an embedding of A into 
MP" 4 followed by a linear projection, may have been introduced into geometric modeling 
in [7J. It is fundamental for Krasauskas's toric Bezier patches [2], and was reworked in the 
tutorial [19]. There, an analog of Theorem 11.111 was formulated. We take this opportunity 
to correct an error in notation. The tautological projection is an algebraic version of the 
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moment map of the toric variety underlying a toric patch, and not the actual moment map 
from symplectic geometry. ^ 

2.3. Rational varieties. We study the algebraic relaxation of our previous notions, re- 
placing the real numbers by the complex numbers so that we may use notions from algebraic 
geometry. A rational patch has a parametrization by rational functions. Clearing their 
denominators gives a new collection {f3 a | a G A} of blending functions which are poly- 
nomials. These polynomials define complex-valued functions on C d , and so the blending 
functions give a map 

fa ■ <£ d \B — > CP* 4 , 
which is defined on the complement of the set B where all the blending functions vanish. 
This is called the base locus of the map (3c- This map (3c extends the map (3: A — > Xp. 
We write (3 for /?c, P" 4 for CP" 4 , and in general use the same notation for maps defined on 
complex algebraic varieties as for their restrictions to subsets of their real points. 

As with linear projections, we write (3 : C d > P" 4 to indicate that (3 is only defined on 

C d \ B. Such a map between algebraic varieties that is given by polynomials and defined 
on the complement of an algebraic subset is called a rational map. 

Let Yp be the Zariski closure of the image of (3. One important consequence of the patch 
being rational is that Xp is a full-dimensional subset of the real points of the algebraic 
variety Yp. In particular, Xp is defined locally in MP- 4 by the vanishing of some polynomials, 
and this remains true for any image F(A) of a map F (11.11) given by our original blending 
functions and any choice of control points. This is important in modeling, for these implicit 
equations are used to compute intersections of patches. 

Suppose that Xp has rational linear precision and that (3 = {(3 a | a G ^4} are rational 
blending functions with linear precision. By Proposition 11.61 the composition 

A ^ Xp -- A ^ A 

is the identity map. Since A has full dimension in IR d , it is Zariski dense in C^, and so the 
composition of rational functions 

(2.8) C d — Yp -^-> C d 

is also the identity map (where it is defined). 

We introduce some terminology to describe this situation. A rational map <p : Y > Z 

between complex algebraic varieties is a birational isomorphism if there is another map 

ip : Z > Y such that tp o if> and ip o tp are the identity maps where they are defined. In 

particular, Y has a Zariski open subset U and Z a Zariski open subset V such that <p\u is 

one-to-one between U and V. The projection tt : C > L of Figure [T] is not a birational 

isomorphism. Indeed, if x G C is a point in 7r _1 (p), then p, E = x, E and this plane meets 
the complement C \ B of the base locus in x and one other point. 

Thus if the composition (12.81) is the identity map, then the map tta is a birational iso- 
morphism from Yp > CP d . We deduce the algebraic-geometric version of Proposition [L6l 

Theorem 2.9. A patch [3 = {f3 a \ a G A} has rational linear precision if and only if the 
complexified tautological projection 

7U : Yp ► CP d 
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is a birational isomorphism. 

Using well-known properties of birational projections, we use this characterization of 
linear precision to deduce necessary conditions for a patch to have linear precision. A 
rational variety of dimension d is one that is birational to C d . 

Theorem 2.10. // a patch (3 = {j3 a | a e A} has rational linear precision, then 

(1) Yp is a rational variety, 

(2) almost all codimension d planes L containing the center of the tautological 
projection meet Yp in at most one point outside of E A . 

By condition (2), Yp has an exceptionally singular position with respect to Ej±. Typi- 
cally, Yp does not meet a given codimension d+1 plane and its intersection with a given 
codimension d plane consists of deg(Yg) points, counted with multiplicity By Condition 
(2), not only does Yp meet E A , but if L has codimension d and contains E^, then most of 
the deg(Yg) points in Yp D L lie in Ej,. We will see this in Examples 12.111 and 13. 161 

Proof. Statement (1) is immediate from Theorem 12.91 

Since the birational isomorphism 7r^ : Yp ► CP d is the restriction of a linear projec- 
tion with center E^, then for almost all points y G Yp \ E^, we have 

(jJT^a) H Yp = {y} U (E A H Yp) . 

Since all codimension d planes L that contain E A have the form (y, EX) for some y ^ E A 
and almost all meet Yp, statement (2) follows. ^ 

Consider the geometric situation of Figure [TJ Suppose now that E is the line tangent to 
C at the point B. Then E fl C is the point B, but counted with multiplicity 2. A linear 
projection with center E restricts to a birational isomorphism of C with P 1 . 

2.4. Unirational varieties. The necessary condition of Theorem 12 .10( 1) gives an impor- 
tant but subtle geometric restriction on patches that have rational linear precision. 

The map (3 provides a parametrization of an open subset of Yp by an open subset of 
C d . Such parametrized algebraic varieties are called unirational. Unirational curves and 
surfaces are also rational, but these two notions differ for varieties of dimension three and 
higher. Clemens and Griffiths [2] showed that a smooth hypersurface of degree 3 in P 4 is 
not rational (these were classically known to be unirational, e.g. by Max Noether). 

Thus parametric patches of dimension at least 3 will in general be unirational and not 
rational. This does not occur for Bezier simploids or toric patches. 

Example 2.11. Consider again the Wachspress coordinates f ll.lOp of Example 11.91 which 
have linear precision. Let Yp C P 5 be the image of the blending functions fll.lOp . We study 
the base locus of the tautological projection 7r^ of this patch. We first find equations for 
Yp as an algebraic subvariety of P 5 . 

If we divide Z^- 1 ]/^ ] + oj + Z^ 1 ]/^ °] by ^ ne product of the linear forms defining 

the hexagon we get 



(1 - x)(l + y - x) + 2(1 + x)(l - x) + (1 + x)(l +x-y) 



= 4 - 2xy , 
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(2.12) fj [zl] ,3^ + + V{ _„ } = ffaffa + 2,3^ + ,3^ 
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If we let (y a \ a E A) be natural coordinates for F A , then (12.121) gives the quadratic 
polynomial which vanishes on Yp, 



2/[ I ^[0]+2|/ [ op [ _o ] +2/ [ ip [ _o ] 



Cyclically permuting the vertices of the hexagon gives two other quadratics that vanish on 
Yp, but these three sum to 0. There is another quadratic polynomial vanishing on Yp, 



There is an additional cubic relation among the blending functions, which gives a cubic 
polynomial that vanishes on Yp 



(2.13) 



These relations, three independent quadratics and one cubic, define Yp as an algebraic 
subset of P" 4 = P 5 . The quadratic equations define Yp, together with the 2-dimensional 
linear space cut out by the three linear equations. 



[2.14) 









y \--\] + ^ + y\u + y[\] + v[jn + v[ 
y [\]- y [--\] + y[i]-y[-i] 



y [\] 



y [--\] + yjs\ 



y [A) 



Thus Yp has degree 7 = 2 3 — 1. A general codimension 2 plane L containing the center 
Ea of the tautological projection tt^ will meet Yp in at most 1 point outside of E^, by 
Theorem 12.10( 2). Since L (1 Yp has degree at least 7, the other 6 points must lie in the 
base locus E^ fl Yp, which is in fact a reducible cubic plane curve. 

The three linear forms (12.141) also define the tautological projection, n^, so this linear 
subspace is the center of projection E^, and the base locus B is defined in E^ by the 
cubic (I2.13P defining Yp. Let us parametrize the center E^ as follows. For [a, b, c] G P 2 , set 



y [-l] 


:= a + b 


y [\] 


:= —a + b 




:= a — c 


y U\ 


:= —a — c 










y \h\ 


:= a — b + c 


y [-l) 


:= —a — b + c 



Then we have 



y [z\] y $} m 



2a(a 2 



b 2 + bc- c 2 



LOJ L1J i-u L o 

the product of a linear and a quadratic form. 
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3. Toric Patches 



Toric patches, which were introduced by Krasauskas [H], are a class of patches based 
on certain special algebraic varieties, called toric varieties. For a basic reference on toric 
varieties as parametrized varieties, see the book [20], particularly Chapters 4 and 13, and 
the articles [U [19], which are introductions aimed at people in geometric modeling. We 
first define toric patches, give some examples, reinterpret Theorem 12.91 for toric patches, 
and then state the main open problem about linear precision for toric patches. 

3.1. Toric Patches. We regard elements a = (ai, . . . , ad) G Z ' as exponent vectors for 
monomials in the variables x — (xi, . . . , xj) 

rr a ■— T a U 2 . . . ™ a d 
ju . — Ju Ju 2 d 

Let A C Z d be a finite set of integer vectors whose convex hull A is a polytope of dimension 
d in R d . This implies that A affinely spans M. d . Suppose that w = {w a G K> | a G A} is a 
set of positive weights indexed by A. These data (A,w) define a map <£a,w' (C x ) d — > P" 4 , 

(3.1) (fA, w ■ x i — ► [w a x a I a G A] ■ 

Let Ya,w be the closure of the image of (C x ) d under the map (fA,w When the weights w a 
are equal, this is the toric variety parametrized by the monomials in A. In general, it is 
the translate of that toric variety by the element w of the positive real torus R>, which 
acts on P- 4 by scaling the coordinates. 

Let X^ w C MP- 4 be the closure y?_4 jU ,(IR>). This is the non-negative part of the translated 
toric variety Y^ w . That is, if IRP> C P" 4 is the set of points having non-negative (real) 
homogeneous coordinates, then X^ w = Y^ w H MP>. It is homeomorphic to the convex 
hull A of the vectors in A, as a manifold with corners [10J § 4]. 

The set A is primitive if it affinely spans Z d . We may assume that A is primitive, for 
there is always some primitive set A' C Z d with Y_^ )W = Y^ )W . First note that translating 
every vector in A by a fixed vector a' multiplies each coordinate of the map ip^,w by the 
monomial x a ' , which does not change tfA,w as a map to projective space. Thus we may 
assume that A contains the zero vector. Since A affinely spans M d , it generates a subgroup 
TLA of Z d which is isomorphic to Z d . Let b 1; . . . , G Z^4 be the basis corresponding to 
the basis of Z d under an isomorphism Z^4 7h d . If we let A' C Z d be the image of A 
under this isomorphism and B := {bi, . . . , b^}, then the map <Pa, w factors 



where \b is the surjective map that sends x to (x bl , . . . ,x hd ). Note that (Pa',w is injective 
but that xb has fibers of size |Z d /Z^4|. Thus <pa,w is not injective unless A is primitive. 

Definition 3.2. A toric patch of shape (A, w) is any patch {/3 a | a G ^4} such that the 
image, X@, of A under the map (12.41) given by the blending functions is equal to the 
non- negative part X^ w of the translated toric variety Ya,w @ 



(C 



Xb 





Remark 3.3. Our geometric perspective that the fundamental object is the image Xp 
requires us to absorb the weights w into our definition and distinguishes patches with 
different choices of weights. /~) 
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Example 3.4. Let i:= {0,l,...,fi}cZ and let w G be any set of weights. Then 
the map (Pa,w is 

X I > [w , WiX, W2X 2 , . . . , w n x n ] . 

The map x 1— > sends A = [0, n) to [0, 00). Precomposing <p^ tW with this map gives a 
toric patch of shape (A, w). 

[0,n] 9 x ^ [w , wi(^), ^ 2 (^) 2 , ^n(^) n 

= [tfo(n — x) n , — x) n , W2X 2 {n — x) n_2 , . . . , w n x n ] . 

Observe that if we choose weights Wi = (™) , substitute ny for x, and remove the common 
factors of n n , then we obtain the blending functions for the Bezier curve of Example II .SI 
Note that replacing x by ny also replaces i G A by ~. ^ 



Example 3.5. Krasauskas [M] generalized the classical Bezier parametrization to any 
polytope with integer vertices. A polytope A is defined by its facet inequalities 

A = {x G W l I hi(x) >0,i = l,...,N}. 

Here, A has N facets and for each i = 1, . . . , N, hi(x) = (vj, x) + q is the linear function 
defining the zth facet, where Vj G Z d is the (inward oriented) primitive vector normal to 
the facet and q G Z. (Compare these to the functions in Example 11.91 ) 

Let A C A fl Z d be any subset of the integer points of A that includes its vertices. Let 
w = {w a I a G A} C M> be a collection of positive weights. For every a G A, Krasauskas 
defined toric Bezier functions 

(3.6) p a (x) := wMx^hix)^ ■■■h N (x) hN{a) . 

For x G A, these parametrize a patch Xp. 

Krasauskas defined a toric patch to be any reparametrization of such a patch. Our 
definition of a toric patch (Definition 13.21) agrees with Krasauskas's, when A = A n Z d 
consists of all the integer points in a polytope: Observe that the map (3 in (13.61) is the 
composition of a map H : A — > M N given by x 1— > (/ii(x), . . . , h^{x)) with a rational map 
V? : ► MP^ given by 

(p : («i, . . . , u N ) 1 — > [WaWj • • ■ Mjv I a G . 

This map <£> factors through the map y?^ )lu . Indeed, define a map 

/ A : (C X ) JV — > (C x ) n 

N 



u = (ui, . . . ,un) 1 — > t = (ti, • • • , in) where tj := u| v * ,e ^ 



=1 



Then 



V V V 

i=l i=l i=l 

And so f(u) = u c ■ <fA,w(fA(u))- This shows that Xp C X^ tW , and Krasauskas [151 §4] 
shows that this inclusion is an equality. /~) 
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Example 3.7. A standard d-simplex of degree n is the simplex in K d , 

nA d := {x G M. d | xi H h x d < n and x { > 0} . 

For a point a = (ai, . . . , a^) in nA^ PI Z d set |a| := ai + • • • + a d and let w a := (™) = 
— ; — 77 — r-m be the multinomial coefficient. This gives a system of weights for A = nAnZ d . 

ai\---ad'.(n— \a\)\ o j o 

Then Krasauskas's toric Bezier patch for nA d has blending functions 

i=l 

If we substitute n?/j for and then remove the common factors of n n , we recover the 
classical Bernstein polynomials [SI §4.1]. ^ 

Example 3.8. The Bezier simploids [6] are toric patches based on products of Bezier 

simplices. Suppose that d±, . . . , d m and ni, . . . , n m are positive integers. Set d :— d\ H h 

d m . Write M. d as a direct sum M. dl © ■ • • © M. dm and suppose that for each i = 1, . . . , m, the 
scaled simplex UiA d . lies in the ith summand M. d '. Consider the Minkowski sum 

A := ni A dl + n 2 A d2 + ■ ■ ■ + n m A dm C R d , 

which is just the product of the simplices niA di . 
If Xj := (Xi s i, . . . , Xi,di) are coordinates on M. di , then 

A = {x = (xi, . . . , x. m ) | x id > and x^ H h x i;di <n { i = 1, . . . ,m} . 

Let A := AnZ d be the set of integer points of A. Given an integer point a = (ai, . . . , a m ) G 
A, set u; a := YllLi (a*)' the product of multinomial coefficients. This gives a system of 
weights, and let Yj^ w be the resulting toric patch. Krasauskas's toric patch of shape A has 
blending functions 

m 

&(x) : = n^w. 

1=1 

where /3 ai (xj) for Xj G M di are the Bernstein polynomials. The resulting patch is called a 
Bezier simploid in [6]. ^ 

We state our main result about linear precision for toric patches, which is a useful 
reformulation of linear precision for toric patches. 

Theorem 3.9. Let A C 7L d be a primitive collection of exponent vectors, w G R> be a 
system of weights, and define the Laurent polynomial 

f = fA,w ■= ^2w a X a . 

a<=A 

A toric patch of shape (A, w) has rational linear precision if and only if the rational function 
ipA,w '■ C ' > C d defined by 

(3-10) T^iTT - /' X2 7r~^> Xd ^~f) 

f V ox i 0x2 ox d / 

is a birational isomorphism. 
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Remark 3.11. The map (I3.10p has an interesting reformulation in terms of toric deriva- 
tives. The zth toric derivative of a Laurent polynomial / is Xi-Jj^-f. The toric differential 
-Dtorus/ is the vector whose components are the toric derivatives of /. Thus the map ipA,w 
is the logarithmic toric differential -Dtoruslog/. 

In the proof below, we show that the map ip A , w comes from a map C d > CP d , 



(3-12) x i — ► / : X!—f : x 2 ^— f : ••• : x d — f] 

L Ox i 0x2 Oxd J 

Let F A>W be the homogenization of / with respect to a new variable xq. Then 

d 

deg(F A>w )F A>w = ^ Xj—F AjW . 

i=o x% 

Thus, after homogenizing and a linear change of coordinates, the map (I3.12p is defined by 
the formula 



D 



torus ^ A,w 



Xqt; — F Aw : x\- — F Aw : ••• : x d - — F Aw ] 
Oxq Ox i Oxd 



This last form of the map (13.121) gives an appealing reformulation of Theorem 13.91 which 
we feel is the most useful for further analysis of linear precision for toric patches. 

Corollary 3.13. The toric patch of shape (A, w) has linear precision if and only if the 
map D^ OTUS F AiW : CP d > CP d is a birational isomorphism. 

Proof of Theorem \3.9i . By Theorem I2.9[ a toric patch of shape (A, w) has linear precision 
if and only if the complexified tautological projection 

*A ■ Y Aw ► CP d 

is a birational isomorphism. Precomposing this with the defining parametrization of Y A>W , 
<PA W '■ (C x ) d <^-> Y AtW (which is a rational map C d > Y A;W ), gives a map 

C d --'^ Y Aw --^ CP d 

( 3 - 14 ) x i — > [w a x a \aeA] i — ► ^w a x a (l,a). 

ae-4 

The initial (Oth) coordinate of this map is 

^ w a x a 

a£A 

which is the polynomial f A>w . For % > the zth coordinate of this composition is 

/ J a i ' Ul a X — X i~Q^~ JA,w ■ 

aeA 1 

If we divide the composite map (13.141) by its Oth coordinate, f AjW , we obtain the map i/) A>w . 
This map is birational if and only if ir A is birational, as <p A)W is birational. ^ 
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Example 3.15. Let A := {0, 1, . . . , n} C Z as in Example EH Let w E M> +1 be any set 
of weights. Then the polynomial / = f^, w is 

f(x) = WQ + WlX^ \-w n x n . 

Suppose that this toric patch has rational linear precision. By Theorem I3.9[ the log- 
arithmic toric differential xf'(x)/f(x) is a birational isomorphism. Thus the fraction 
xf'(x)/f(x) reduces to a quotient of linear polynomials and so f(x) is necessarily a pure 
power, say a(x + a) n with a^O. 

We conclude that the weights are Wi = which are not quite the weights in the 

Bezier curve [L~8l If we rescale x (which does not change Y^ w ), setting y = x/a, then the 
weights become |a| n ("). Removing the common factor of |a| ra , reveals that this toric patch 
is the classical Bezier curve of Example 11.81 ^ 

We close with two problems. 

Problem 1. Classify the toric patches of dimension d that have linear precision. 

The analysis of this section was recently used to classify toric surface patches (d = 2) 
having linear precision [11] . It remains an open problem to understand how to tune a toric 
patch (moving the points A) to achieve linear precision. 

Problem 2. Classify the toric patches of dimension d that may be tuned to have linear 
precision. 

This is a strictly larger class of patches. For example, the classification of [11] shows if 
a toric surface patch has linear precision, then it is triangular or quadrilateral, but it is 
possible to tune a pentagonal patch to achieve linear precision. 

Example 3.16. We discuss Krasauskas's version [15] of Karciauskas's pentagonal M- 
patch [12], which has linear precision. Let A be the integer points in the pentagon Qj 
shown below, and let {w a | a E A} be the indicated weights. 




• • • 3 5 2 

(0,0) (1,0) (2,0) 

Lattice points A weights {w a \ a E A} 

These define a toric patch X^ w and translated toric variety Y^ w . This subvariety of P 
has degree 7 and is defined by 14 quadratic equations, all of the form 

UaUh — VcVd where a, b, c, d E A with a + b = c + d . 

This has a monomial parametrization 

<PA,w ■ (s,t) .— > [s a H a2 | a= (a u a 2 ) E A] . 
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Its composition with the tautological projection tt^ gives a map g : P 2 > P 2 of degree 4, 

with 3 base points at (—1, —1), (—1, — |), and (— |, —1) 

We tune this patch (or rather the projection ttj\) by moving the non-extreme points of 
A within the pentagon. Let B consist of the vertices of the pentagon, together with the 
three non-vertices 



(0,2) 



(04: 



;i,2) 



'8 8^ 
>7' 7 ' 



►(2,1) 



(0,0) (|,0) (2,0) 
Lattice points B 

The composition of <pa,w with this new tautological projection 7i& gives the rational map 

2s{2s + t + 3) 2t{s + 2t + 3) 



> + l)(2s + 2t + 3) ' (t + l)(2s + 2t + 3) 

-,-1), (0,- 



This is given by quadrics and has 3 base points (—1, —1), (0, — |), and (— ~, 0), and so it 
has degree 1 = 2-2 — 3, and is therefore a birational isomorphism. By Theorem 12. 9[ the 
tuned toric patch (with the lattice points A replaced by B) has linear precision. 

Thus there is a rational parametrization f3 : Q — > X^ w whose composition with the 
tautological projection ttb is the identity map. Indeed, Krasauskas [15J gives the following 
modification of the toric Bezier functions and shows that the map TTg o f3 is the identity. 



^1 



(3.17) 



= 3(2- 


s) 2 (3-s-t) 3 (2-t) 2 




= 5(2- 


s)(3-s-t) 2 (2-t) 2 s(3-s- 


2> 


= 5(3- 


|-t)t(2- S ) 2 (3-s-t) 2 (2- 


t) 


= 7(3- 


f - t)t(2 - t){3 - s - t){2 - s] 


s(3 


= 2(3- 


s-t)(2-t) 2 s 2 (3-s- t 1 ) 2 




= 2(3- 


l-t) 2 t 2 (2-s) 2 s(3-s-t) 




= 2(3- 


l-t)t(2-t)s 2 (3-s-l) 2 




= 2(3- 


f - t) 2 t 2 (2 - s)s(3 - s - |) 





This is a modification of the toric Bezier functions for the pentagon. If, in the defini- 
tion (I3.6p . we replace the form s by s(3 — s — |) and the form t by (3 — | — t)t, then we 
get these modified blending functions (13.171) . 
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Let us consider the geometry of linear precision for this tuned patch in the spirit of 
Section [2J The projection ir^ is defined by the three linear forms. 



X 
Y 
Z 



gyg, + 2yg| + ,y {l] + 2y [?] + yg, 
§y [?] + 2y [ o ] + fy E] + 2y^ + y [?] 



Its center £g is defined by the vanishing of these three forms. A general linear subspace L 
of codimension 2 containing Eg is defined by equations of the form 



X 



xZ 



and 



Y 



yZ. 



This subspace meets Y^ w in four points. One point lies outside of E B while three points 
lie in the center Eq. We give the three points. 



2 — -2. 



-5 







— 







-5 7 -2 



multiplicity 1 multiplicity 1 multiplicity 4 

We write the coordinate y a of a point y G MP- 4 at the position of the point a. 



4. Iterative proportional fitting for toric patches 

The toric patch X^ w from geometric modeling appears naturally in algebraic statistics in 
the form of a toric model, which leads to a dictionary between the subjects. We show that 
a toric patch has rational linear precision if and only if it has maximum likelihood degree 
1 as a statistical model. As a consequence, we present a new family of toric patches with 
rational linear precision. Finally, we propose the iterative proportional fitting algorithm 
from statistics as a tool to compute the unique reparametrization of a toric patch having 
linear precision. 

In algebraic statistics, the image of R> under the map (13.11) is known as a toric 

model [T5], §1-2]- They are more commonly called log-linear models, as the logarithms of 
the coordinates of <p^, jW are linear functions in the logarithms of the coordinates of M>, or 
discrete exponential families, as the coordinates of (Pa,w are exponentials of the coordinates 
of R d , which are themselves logarithms of the coordinates of Mz. 

We identify the non-negative orthant MP> with the probability simplex 

A" 4 := {yeM^ | |y| 1}, 
and so we may also regard X^ tW as a subvariety of A A . 
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The tautological map (12.51) appears in statistics. Given (normalized) data q G A" 4 , 
the problem of maximum likelihood estimation is to find the toric model parameters (a 
point in IR>) that best explain the data q. By Lemma 4 in [5], the maximum likelihood 
estimate for the log-linear model X^ w is the unique point p G X^ w such that 

Ta(p) = n A {q) so that p = ^(^(g)). 

Thus inverting the tautological projection is necessary for maximum likelihood estimation. 

Catanese, et. al. pQ defined the maximum likelihood degree of the model Xj^ w to be 
the degree of tt^ 1 , as an algebraic function. Equivalently, this is the number of complex 
solutions to the critical equations of the likelihood function, which is the degree of the 
tautological map tt a from the Zariski closure Y^ w of X^ w to CF d . By definition, a toric 
patch has rational linear precision if n^ 1 is a rational function (an algebraic function of 
degree 1). 

Proposition 4.1. A toric patch X^ w has rational linear precision if and only if X^ w has 
maximum likelihood degree 1. 

An important family of toric models called decomposable graphical models are known to 
have maximum likelihood degree 1 [Tol p. 91]. Therefore, these models, which are not in 
general Bezier simploids, have linear precision. They typically have large d > 2 dimension. 

Darroch and Ratcliff [5] introduced the numerical algorithm of iterative proportional fit- 
ting, also known as generalized iterative scaling, for computing the inverse of the tautologi- 
cal projection n^ 1 . Their algorithm requires that the data A be in a normal, homogeneous 
form. 

Observe first that the toric patch X^ w does not change if we translate all elements of 
A by a fixed vector b, (a i— > a + b) , so we may assume that A lies in the positive orthant 
IR>. Scaling the exponent vectors in A by a fixed positive scalar t G R> also does not 
change X A w as i ^ i ! is a homeomorphism of M> that extends to a homeomorphism of 
R>. Thus we may assume that A lies in the standard simplex A^ in R d : 

A d = {x G R d | x > and \x\ < 1} . 

Lastly, we lift this to the probability simplex A J G 

A+ := {y G R 1+d \yi>0 and \y\ = 1} , 

by 

A3 sl i — > a + := (l-|a|, a) G A+ . 
Since for t G K> and x G M>, 

{t,tx) a+ = t 1 ~ |a| (tx) a = tx a , 

we see that replacing A by its lifted version A + also does not change X^ w . We remark 
that this is just a way to homogenize the data for our problem. 

Since these are afline transformations, the property of affine invariance for patches shows 
that the affine map that is the composition of these transformations intertwines the original 
tautological projection with the new one. 

We describe the algorithm of iterative proportional fitting, which is Theorem 1 in [5] . 
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Proposition 4.2. Suppose that A C and y G conv(A). Then the sequence 

{p(») | 71 = 0,1,2...} 



whose a- coordinates are defined by p^ ■= w a and, for n > 0, 

„(n+l) . = (n) . 



converges to the unique point p G suc/i that vr^(p) = n^(y). 

We remark that if A is not homogenized then to compute ir^ (y) for y G conv(^4), 
we first put A into homogeneous form A + using an affine map ip, and then use iterative 
proportional fitting to compute ir~^ + (ip(y)) = 7T^ (y). We also call this modification of 
the algorithm of Proposition 14.21 iterative proportional fitting. Thus iterative proportional 
fitting computes the inverse image of the tautological projection, which, by Theorem 12.61 
gives the unique parametrization of X^ w having linear precision. 

Corollary 4.3. Iterative proportional fitting computes the unique parametrization of a 
toric patch having linear precision. 

Little is known about the covergence of iterative proportional fitting. We think that it 
is an interesting question to investigate this convergence, particularly how it may relate to 
maximum likelihood degree. 
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